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GENERAL CONSIDERATIONS 


The knowledge of the density distribution functions in the adsorbate 
phase is of extremely big importance, both in theoretical, and in practical 
investigations of physical adsorption. The theoretical calculations are the 
only possibility to obtain this function at the present time, since the 
problem of its experimental determination still remains unresolved. 
Because so far no experimental method exists to solve its evaluation. 

Recently Kuni [1, 2] has developed the elegant formula for surface 
distribution functions, using the methods of functional analysis. The Kuni’s 
method would be a very useful one, if it were not limited by the condit- 
ion, that the method is not applicable for the most interesting region 
near the solid surface. In this paper we extend the Kuni's theory for the 
region near the solid surface in the case of liquid adsorbate phases, for- 
med by simple liquids. 

Following Kuni we consider the liquid adsorbate phase as a liquid in 
the presence of an external field due to the solid surface. The basic Kuni’s 
idea is to consider the adsorbate phase; first in the presence of this exter- 
nal solid surface field e (r ) (assumed to be of one-particle type) and next 
in the presence of an “auxiliary” external field e, ( r). 

Let g™ (m, ..., Tn) and g™ (Ty, «nes Tn) are the n-body Mayer-Ursell's cor- 
relation functions in the adsorbate phase, being in the presence of the 
external solid field e (r), and in the presence of the "auxiliary” field 


Eo (r), respectively. Then the Kuni's method leads to the following ex- 
pansion: 
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where: V is the volume of the adsorbate phase; next Ae (r)=e(r)— 
toe | r); k and T are the Boltzmann constant and absolute temperature, 
respectively. The above expansion giving the Mayer-Ursell's correlation 
functions is simultaneously a a kind of expasion for common symmetrical 
distribution functions gt” (r;, Bs. Tn) because of the simple connections 
between themselves. For example: 


gT) = pM (7) (2) 
IPE, Te) = PATE Te PM) PM) (3) 


Later, in the case of 9 (1) in a liquid adsorbate phase one may neglect 
in the Kuni's expansion all terms with n+m > 2, since three and higher 
body correlations in liquids contribute in a little measure to the total 
properties of liquids [3]. 

Thus, in our case of the density distribution function ọ%® (ri) the Ku- 
ni’s expansion reduces to the following formula: 


p"imy= [exp AEE Je" j+ | [aera p re pe] « 
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where 0, (r, bens Th) are the common symmetrical distribution functions 
in the adsorbate phase in the presence of the auxiliary” external field 
Eo (r). According to the Kuni's method the ’auxiliary” external field 
may be chosen quite arbitrary, but for the practical reason it should be 
chosen so that the auxiliary correlation functions may be easily evalua- 
ted. 

In an original procedure developed later by Kuni, it is possible only 
for the region of the adsorption space, where the reciprocal of the distance 
from the surface can be considered as a small parameter. In this way only 
the asymptotical form of these distribution functions has been obtained 
by Kuni and co-workers. In this paper we are going to propose a new 
method for calculating the auxiliary field and the auxiliary distribution 
functions. The method is based on the tunnel theory of dense fluids [4, 
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5], and is applicable for the whole adsorption space, in particular for the 
region near the solid surface. 


The tunnel theory follows the basic idea of the cell theory that one 
should divide the whole liquid into subsystems, sufficiently simple to per- 
mit calculation of their properties. Thus, the liquid volume is divided into 
subsystems being hexagonal cylinders (or tunnels), arranged in close- 
-packing. Next, one assumes that each tunnel contains an equal number of 
molecules. Since each tunnel contains a large number of molecules, it is 
not necessary to take into account the fluctuations of the number of mole- 
cules in each tunnel. The molecules in different tunnels move independen- 
tly, but only in their own tunnels, the walls of each tunnel being for- 
med by molecules of the neighbouring tunnels. 


In using the tunnel model to our purpose, we shall approximate the 
solid surface by a plane (it will be the plane) and assume the adsorption 
potential to depend on the distance from the surface (xy plane), only. 
The above two assumptions are usually made in the majority of public- 
ations done in this field so far. 


Now consider a hypothetical auxiliary liquid phase, consisting of the 
same kind of molecules, and being enclosed in the same volume V as our 
liquid adsorbate phase. Let the density of the auxiliary phase be equal to 
the density @ of our adsorbate phase at long distances from the solid sur- 
face (in the bulk liquid). We divide the auxiliary phase into the hexagonal 
tunnels so that the axes of these tunnels should be perpendicular to the 
plan xy bounding partially the auxiliary phase. Next we suppose the 
existence of an auxiliary external field £o (T) acting on the auxiliary 
phase so that in spite of symmetry lack of interactions near the xy plane, 
the density should be still equal to ọ through the auxiliary phase. Calcu- 
late the auxiliary external field, to use it in the eqn (4). 


Following Barker we consider the molecules of neighbouring tunnels 

s ”smeared” uniformly on cylindrical surfaces, whose radiuses and sur- 
iis densities are tabulated in Barker’s papers [2, 3]. Let us calculate the 
interaction energy V (m) of molecule at the point rę with the _walls of 
these cylindrical surfaces. To this purpose we transform vector m to cy- 
lindrical coordinates T (i, z Following Barker we get: 


A L 2m 
VF) = VI; zi) = Y faz jag Ry ds Uld) (5) 
sai 0 o 
where: L is the length of the tunnel; o, is the surface density in the s-th 
cylindrical surface, defined as the ratio of the number of molecules in the 
s-th cylindrical surface, to this surface area; Rs is the radius of the s-th 
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cylindrical surface. Next: dż = (z, — 2;)?- R2,-- r2,—2r, Rs cos Q; and U (d) is 
the interaction energy between two molecules at the distance d. 

Let us denote by Voo the value of V for the distance from the xy 
plane, where the symmetry lack of interactions may be neglected. Thus, 
the potential V may be considered as consisting of two parts; Voo and 
AVe, the second of which being due to the symmetry lack of interactions 
near the xy plane. 


BIZ) = Wri) + AVE (Ziri) (6) 


However, the symmetry lack of interactions near the boundary xy ari- 
ses not only from the symmetry lack of interactions with the cylidrical 
surfaces, but from the symmetry lack too, in the interactions between 
the molecules being in the same tunnel. Let us denote the average inter- 
actions energy of a considered molecule with all other molecules in the 
same tunnel by AVr. It is done by the following obvious equation: 


AV,(7) » (UA Bg a. Bay = 


Above the interaction space Vr is simply the volume of the tunnel. 

The function Eo! (Te, r, T3) has the following meaning. If there is a mo- 
lecule at the point r,, then the density of molecules at the point r in the 
same tunnel is equal to gl ( 4, T) T;). For our auxiliary phase having the 
density o through the whole space, the following simple relation can be 
written: 


ROR) = PĘD) ©) 


Evaluating AV, from eqn. (6) and AVr from eqn. (7), we get the 
required expression for the auxiliary external field. 


KIT = Gizi = — [Ar (zur) + AVIZA)] (9) 


The above expression | can be inserted into eqn. (4). Now, it is seen that 
in addition to s (m) the two-body distribution function go?! (7, Tj) 
must be known, to use eqn. (4) for our purpose. However, the method of 
evaluating go?! ( Ty Tj) in the tunnel space has been in detail described in 
Barker’s paper [6]. 

Now, we have all the functions whose knowledge is required when 
using eqn. (4). With the assumptions concerning the auxiliary phase the 
eqn. (4) may be rewritten to the following simpler form: 


pz) = aa aem pe[ [eran (4642) -] on} (10) 
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As it should be expected, the density KUCH ) in the adsorbed phase is 
a function of the adsorbate density o in the bulk phase. 


APPROXIMATE MODEL 


In this section we shall apply the eqn. (10) to hard-sphere fluid, 
being in the presence of an external field which is defined as follows: 


+œ for 0 pa z<a 
E(FF)=€lz)=<Emia tor aCz<b (11) 
0 for bKz<+wo 


It is well-known that the hard-sphere fluid is a good approximation for 
many liquids having the slow-collision diameter equal to the particle 
diameter of this hard-sphere fluid. Next, it is known that the function in 
eqn. (11) with suitable chosen parameters; a, b, tmin is a good approxi- 
mation for majority of adsorption potentials. In addition, the application 
of hard-core interactions in theoretical descriptions of the properties of 
fluids, leads to extremely simple results. For this reason we expect to 
get simple formulas for 9“ ( r; ) in our adsorbate phase, when considering 
it as a hard-sphere fluid in the presence of the external field from eqn. 
(11). Next, we shall later assume: 


+oo for OKd<D 
U(d) = (12) 
0 for DKd<two 


Above D is the slow-collision diameter. 

In the case of the hard-sphere fluid, the interactions of a considered 
molecule with all cylindrical surfaces are equal to zero, except for the 
interactions with the first nearest cylindrical surface. This interaction 
is described as follows: 


(13) 


— 


0 for OK nN<IR-D) 
VIz, ) = 
+œ for (R-DJK rn +œ 


Here R is the radius of the first cylindrical surface, its value being 
connected with the density ọ of the adsorbate [4]. 


4 
A, Sale (14) 
r= (r) 
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For really existing fluids (liquids) D is about 90% of R. (R=1, 11D for 
liquid argon in 84°K [6]). It justifies our next approximation that, 


+œ for 0<¢1Z-ZI1<D 
viai=unzezi= | 53 (15) 
0 for DIZ-Zj|C+0o 
Next, we approximate the function g.?! (TT) by the following one: 
0 for OKIz-ziI D 
p tor D<Izi-zji <+ 


For a variety of practical situations D is much higher than (b—a). In 
the case of gas-solid adsorption system of carbon tetrachloride adsorbed 
en graphit, the value (b—a) found by us was about 107*cm. We shall 
therefore assume (b—a) < D. Next, for majority adsorption systems D>a. 
For example D=3,4 A, and a=2,9 A for argon on carbon black [7, 8]. 
Thus, we shall later assume D > a. With these all above assumptions eqn. 
(4) yields: 


tafor OC Z<a 

p"F)=0 
2.tor agz2<b 

= = | exp je | [Pe aA,-p *(b-a) Ay ©xP Enis | 
3.for bKz Ka+D : 


p(T 7) = pM z,2) = | (16) 


na. p+p a-z) Ar — pib- -g) A; exp Emin 


4. for a+Dz<b+D 


(17) 


Emin 


pr) = p+p (D-z+b] Ar exp ÍRT 


5 for b+DKZK+w 
pF) =p 


Above Ar is the free cross-sectional area of the tunnel being equal to 
x (R—D)?. 


Let us denote by N the number of adsorbed molecules (the nett Gib- 
bses adsorption). Then: 


N= len - pjar = Be | trip + eme] (18) 


Above S is the adsorbent area, and 
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Emin 
AT)=Arlb-a) [exl £) - i (19) 
In the above equation we shall approximate the value, since for the 
[exp] -1] by exp See 


majority of adsorption systems the last value is about 10° at room 
temperatures. Next: 


(20) 
BIT) =alb- o]Ar|exp nie Ere |--(b-afa, [exp 26 | + ala+0-b)A, - 


-|-o)la+D-BJA,|exp Sz — [301° - $A HOrolb-a1,|exp Eni- 
ma ieran 


From the experimental plot ME "against go one can get both a, and 


B. Let us insert eqn. (19) into eqn. (20), found in the temperature Tn. 


Th) _ (b+D) , (a+D b+Df_ la+D. 
PIL) =A) [2 4l- Ibo „| iD) |+aio+D-b1a,- [230 le50] Ja. 


irs (21) 
By plottin N Ae against ọ in a number of temperatures one gets the 
y 8 2S 


same number of different equations of type (21). From these equations 
one can evaluate the parameters: a, b, Ein, D, S, by means of the method 
of the smallest areas. 

For many liquids whose molecules have not spherical symmetry, the 
hard-sphere model may be used by introducing the concept of "effective 
slow-collision diameter”, and "effective adsorption potential”. These are 
the appropriate potentials, averaged over all coordinates describing the 
state of molecule, except for the center of mass. Here by comparing the 
effective slow-collision diameter in the adsorbed, and in the bulk phase, 
one may extract additional interesting informations about the adsorption 
model. 
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STRESZCZENIE 


W pracy dokonano próby rozszerzenia teorii Kuniego funkcji dystry- 
bucji w fazach powierzchniowych celem opisania obszaru leżącego blisko 
powierzchni ciała stałego. Z tego względu zaaplikowano pewne idee z te- 
orii tunelowej Barkera opisującej gęste płyny. 


PE3IOME 


Nytem paciumupeHua Teopuu noBepxHOCTHbiX (byHkuiMA pacnpeąeneHna 
Kynn, aBTOpbi NbiTałoTca OnucaTb pacnpeĄeneHue NNOTHOCTM BÓNK3u noBepx- 
HOCTM TBepĄoro Tena. Ana 3Toh yenn npHMEHA4IOTCA HEKOTOpPbie ngen TEOPMH 
Bapkepa, onncbiBatoluen NNOTHbIE XKMĄKOCTM. 


